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ABSTRACT 

Context. Orbital obliquity is thought to be a fundamental parameter in tracing the physical mechanisms that cause the 
migration of giant planets from the snow line down to roughly 10 -2 au from their host stars. We are carrying out a 
large programme to estimate the spin-orbit alignment of a sample of transiting planetary systems to study what the 
possible configurations of orbital obliquity are and wether they correlate with other stellar or planetary properties. 
Aims. We determine the true and the projected obliquity of HAT-P-36 and WASP-ll/HAT-P-10 systems, respectively, 
which are both composed of a relatively cool star (with effective temperature T e ff < 6100 K) and a hot-Jupiter planet. 
Methods. Thanks to the high-resolution spectrograph HARPS-N, we observed the Rossiter-McLaughlin effect for both 
the systems by acquiring precise (3 — 8ms~ 4 ) radial-velocity measurements during planetary transit events. We also 
present photometric observations comprising six light curves covering five transit events, which were obtained using three 
medium-class telescopes. One transit of WASP-ll/HAT-P-10 was followed contemporaneously from two observatories. 
The three transit light curves of HAT-P-36 b show anomalies attributable to starspot complexes on the surface of the 
parent star, in agreement with the analysis of its spectra that indicate a moderate activity (log_R(j K = —4.65dex). By 
analysing the complete HATNet data set of HAT-P-36, we estimated the stellar rotation period by detecting a periodic 
photometric modulation in the light curve caused by star spots, obtaining P TO t = 15.3 ± 0.4 days, which implies that 
the inclination of the stellar rotational axis with respect to the line of sight is i* = 65° ± 34°. 

Results. We used the new spectroscopic and photometric data to revise the main physical parameters and measure the 
sky-projected misalignment angle of the two systems. We found A = —14° ± 18° for HAT-P-36 and A = 7° ± 5° for 
WASP-ll/HAT-P-10, indicating in both cases a good spin-orbit alignment. In the case of HAT-P-36, we also measured 
its real obliquity, which turned out to be -ip = 25jj 38 degrees. 

Key words, stars: planetary systems - stars: fundamental parameters - stars: individual: HAT-P-36 - stars: individual: 
WASP-ll/HAT-P-10 - techniques: radial velocities - techniques: photometric 
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1. Introduction 


Since its discovery dMavor &; Quelozl Il995h . hot Jupiters 
have challenged the community of theoretical astrophysi¬ 
cists to explain their existence. They are a population of 
gaseous giant extrasolar planets, similar to Jupiter, but re¬ 
volving very close to their parent stars (between ~ 0.01 and 
0.05 au), which causes them to have orbital periods of few 
days and high equilibrium temperatures (between ~ 400 
and 2750 K]3- According to the generally accepted theory 
of p la net formationlwe refer the reader to the review of 
iMordasini et al.l [201411 . hot Jupiters are thought to form 
far from their parent stars, beyond the so-called snow line , 
and then migrate inwards to their current positions at a 
later time. However, the astrophysical mechanism that reg¬ 
ulates this migration process, as that which causes planets 
moving on misaligned orbits, is still a matter of debate. 

Spin-orbit obliquity, ip, that is the angle between a 
planet’s orbital axis and its host star’s spin, could be the 
primary tool to understand the ph ysics b ehind the migra¬ 
tion process of giant planets (e.g., iDawso iil 120141 and ref¬ 
erences therein). A broad distribution of ip is expected for 
giant planets, whose orbits are originally large (beyond the 
snow line) and (being perturbed by other bodies) highly ec¬ 
centric, and then become smaller and circularize by planet’s 
tidal dissipation during close periastron passages. Instead, 
a much smoother migration, like the one through a proto¬ 
planetary disc, should preserve a good spin-orbit alignment 
and, assuming that the axis of the proto-planetary disk re¬ 
mains constant with time, imply a final hot-Jupiter popu¬ 
lation with ip cz 0. 

The experimental measurement of ip for a large sam¬ 
ple of hot Jupiters is mandatory to guide our comprehen¬ 
sion of such planetary evolution processes in the mean¬ 
ders of theoretical speculations. Unfortunately, this quan¬ 
tity is very difficu lt to determine, and only a few measure¬ 
ments exist (e.g. , iBrothwell et al.l[2014 iLendl et~aT1l20Ui 
iLund et aPl2014fl . However, its sky projection, A, is a quan¬ 
tity which is commonly measured for stars hosting tran¬ 
siting exoplanets through the observation of the Rossiter- 
McLaughlin (RM) effect with high-precision radial-velocity 
(RV) instruments. Such measurements have revealed a wide 
range of obliquities in both early- and lat e-type stars (e.g., 
iTriaud et al.l l201d : lAlbrecht et alJ 20121J). including pl an- 
ets revolving p erpendicular (e.g.. lAlbrec ht, et al.l l2012al) or 
retrogade (e.g. . Anderson et al. [20 1 Ot iHebrard et al* 1 201 it 
lEsposito et al .1120141) to the direction of the rotation of their 
parent stars. The results collected so far are not enough to 
provide robust statistics and more investigations are needed 
to infer a comprehensive picture for identifying what the 
most plausible hot-Jupiter migration channels are. 

Within the framework of GAPS (Global Architecture 
of Planetary Systems), a manifold long-term observa¬ 
tional programme, using the high-resolution spectro¬ 
graph HARPS-N at the 3.58 m TNG telescope for sev¬ 


* Based on observations made with ( i ) the Italian 3.58 m 
Telescopio Nazionale Galileo at the Observatory of Roque de los 
Muchachos, (ii) the Cassini 1.52 m telescope at the Astronomical 
Observatory of Bologna, (in) the Zeiss 1.23 m telescope at the 
Observatory of Calar Alto, and the IAC 80 cm telescope at the 
Teide Observatory. 

1 Data taken from TEPcat (Transiting Extrasolar Planet 
Catalogue), avail able at www.astro.keele.ac.uk/jkt/tepcat/ 
(ISouthworthl 12111 1). 


eral s emesters (iDesidera et alJ 2013. 2 0141 : iDamasso et al.l 
1201,51) . we are carrying out a subprogramme aimed at study¬ 
ing the spin-orbit alignment of a large sample _of known 
trans iting extrasolar planet (TEP) systems dCovino et al.l 
120131) . This subprogramme is supported by photometric 
follow-up observations of planetary-transit events of the 
targets in the sample list by using an array of medium- 
class telescopes. Monitoring new transits is useful for get¬ 
ting additional information, as stellar activity, and for con¬ 
straining the whole set of physical parameters of the TEP 
systems better. When possible, photometric observations 
are simultaneously performed with the measurement of the 
RM effec t, benefitin g from the two-telescope observational 
strategy (|Ciceri et al.l 1 20131) . 

Here, we present the results for two TEP systems: HAT- 
P-36 and WASP-ll/HAT-P-10. 

The paper is organized as follows. In Sect. ll.il we briefly 
describe the two systems, subjects of this research study. 
The observations and reduction procedures are treated in 
Sect.[2j while Sect.[3]is dedicated to the analysis of the pho¬ 
tometric data and the refining of the orbital ephemerides; 
anomalies detected in the HAT-P-36 light curves are also 
discussed. In Sect.|J] the HARPS-N time-series data are 
used to determine the stellar atmospheric properties and 
measure the spin-orbit relative orientation of both the sys¬ 
tems. Sect. [5] is devoted to the revision of the main physi¬ 
cal parameters of the two planetary systems, based on the 
data previously presented. General empirical correlations 
between the orbital obliquity and various properties of TEP 
systems are discussed in Sect. [6] The results of this work are 
finally summarised in Sect. [7] 

1.1. Case history 

HAT-P-36, discovered by the HAT-Net survey ( Bakos et al.l 
[20ll . is composed of a V = 12.5 mag, G5V star (this 
work), similar to our Sun, and a hot Jupiter (mass ~ 
1.8Mj up and radius ~ 1.3i?j up ), which revolves around its 
parent star in nearly 1.3 days. The only follow-up study 
of this system has provided a refinement of the tran¬ 
sit ephemeris based on the photometric observation of 
two planetary transits p erformed with a 0.6 m telescope 
dMacieiewski et ;dl 1 201 31) . 

The discovery of the WASP-ll/HAT-P-10 TEP sys¬ 
tem was almost c oncurrently annou nced by the WASP 
(IWest et all I2009P and HAT-Net dBakos et alJ [2009) 
teams. It is com posed of a V = 11.9 mag, K3V star 
(lEhrenreich fe Desertl 1201 3 ) (mass ~ 0.8 Mq and radius 
~ 0.8 Rq) and a Jovian planet (mass ~ 0.5 Mj up and ra¬ 
dius ~ 1 i?j up ), moving on a 3.7 day circular orbit. Four ad¬ 
ditional tra nsit ob servations of this target were obtained by 
ISada et ahl ( |2012i ) in z * l and J bands with a 0.5 nr and 2.1m 
telescope, respectively, which were used to refine the tran¬ 
sit ephemeris and phot ometr ic parameters. In a more re¬ 
cent study (IWang et al.ll2014) . an additional four new light 
curves, which were obtained with a 1 nr telescope, were pre¬ 
sented and most of the physical parameters were revised, 
confirming the results of the two discovery papers. 


2. Observation and data reduction 

In this section we present new spectroscopic and photo¬ 
metric follow-up observations of HAT-P-36 and WASP- 
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Fig. 1. Phase-folded RV data of a partial transit of HAT- 
P-36b observed with HARPS-N. Superimposed is the best¬ 
fitting RV-curve model. The corresponding residuals are 
plotted in the lower panel. Phase 0 corresponds to the time 
of the planet passing the periastron. 



Orbital phase 


Fig. 2. Phase-folded RV data of a complete transit of 
WASP-11 b observed with HARPS-N. Superimposed is the 
best-fitting RV-curve model. The corresponding residuals 
are plotted in the lower panel. 


ll/HAT-P-10. For both the planetary systems, the RM ef¬ 
fect was measured for the first time and three complete 
transit light curves were acquired. 

2.1. HARPS-N spectroscopic observations 

Spectroscopic observations of the two targets were per¬ 
formed using the HARPS-N (High Accu racy Radial veloc¬ 
ity Planet Searcher-North: iCosentino et al]l2012lf spectro¬ 
graph at the 3.58 m TNG, in the framework of the above- 
mentioned GAPS programme. 

A sequence of spectra of HAT-P-36 was acquired on the 
night 2013 February 21 during a transit event of its planet, 
with an exposure time of 900 sec. Unfortunately, because 
of high humidity, the observations were stopped before the 
transit was over and thus only nine measurements were 
obtained that cover a bit more than the first half of the 
transit. 


Table 1. HARPS-N RV measurements of HAT-P-36 and 
WASP-ll/HAT-P-10. 


BJD(TDB) 

RVfms- 1 ) 

Error (m s 1 ) 

HAT-P-36: 

2 456 345.583549 

-16200.7 

5.5 

2 456 345.594369 

-16213.9 

5.8 

2 456 345.605189 

-16227.2 

4.6 

2 456 345.616014 

-16249.9 

4.6 

2 456 345.626843 

-16240.2 

5.5 

2 456 345.637659 

-16247.8 

5.3 

2 456 345.648483 

-16266.6 

7.4 

2 456 345.659304 

-16306.0 

5.5 

2 456 345.670130 

-16352.7 

6.4 

WASP-ll/HAT-P-10: 


2 456 933.502525 

4907.4 

8.0 

2 456 933.509760 

4901.6 

6.5 

2 456 933.516999 

4907.8 

4.5 

2 456 933.524234 

4902.7 

3.8 

2 456 933.531469 

4905.5 

3.7 

2 456 933.538699 

4910.1 

3.5 

2 456 933.545929 

4902.3 

4.1 

2 456 933.553165 

4907.2 

3.5 

2 456 933.560400 

4908.2 

4.0 

2 456 933.567644 

4902.8 

3.5 

2 456 933.574884 

4913.1 

3.3 

2 456 933.582122 

4913.4 

3.9 

2 456 933.589366 

4907.7 

3.7 

2 456 933.596598 

4903.3 

4.2 

2 456 933.603834 

4904.6 

4.0 

2 456 933.611071 

4899.7 

4.2 

2 456 933.618316 

4891.2 

4.2 

2 456 933.625565 

4887.1 

3.3 

2 456 933.632809 

4880.9 

3.5 

2 456 933.640053 

4879.5 

3.5 

2 456 933.647298 

4877.8 

3.2 

2 456 933.654533 

4877.9 

3.8 

2 456 933.661763 

4879.4 

3.7 

2 456 933.669002 

4885.9 

3.7 

2 456 933.676238 

4887.0 

3.8 

2 456 933.683469 

4891.0 

4.1 

2 456 933.690700 

4884.7 

3.7 

2 456 933.697935 

4881.6 

3.8 

2 456 933.705166 

4884.5 

3.5 

2 456 933.712397 

4873.5 

4.0 

2 456 933.719632 

4881.2 

4.3 

2 456 933.726867 

4879.8 

4.2 


A complete spectroscopic transit of WASP-ll/HAT-P- 
10 b was recorded on 2014 October 2. A sequence of 32 
spectra was obtained using an exposure time of 600 sec. 

The spectra were reduced using the latest version of the 
HARPS-N instrument Data Reduction Software pipeline. 
The pipeline also provides rebinned ID spectra that we used 
for stellar characterization (see summary Table[5] and [9]), 
in particular to estimate the stellar effective temperature, 
T e g, and metal abundance [Fe/H] (see Sect. 14. 111 . Radial 
velocities were measured by applying th e weighted cross¬ 
correlation funct ion (CCF) method ( Baranne et al.1 If 99 til 
iPene et aT1l2002t) and using a G2 and a K5 mask for HAT- 
P-36 and WASP-ll/HAT-P-10, respectively. They are re¬ 
ported in Tableland plotted in Figs.[l]and[2j in which we 
can immediately note that the RM effect was successfully 
observed. It indicates low orbital obliquity for both the sys¬ 
tems. 
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Fig. 3. Light curves of HAT-P-36b transits compared with the best PRISM+GEMC fits. The dates, the telescopes and the 
filters used for each transit event are indicated. Residuals from the fits are displayed in the right panel. Starting from the 
top, the first three light curves (this work) present anomalies that we interpret as occultations of starspot complexes by 
the planet. The third ligh t curve is rela ted to the transit event successive to that of the second one. The last light curve 
was extracted from iBakos et ah (j2012f ) and reported here for comparison. Here phase 0 corresponds to the mid-time of 
the transit. 


2.2. Photometric follow-up observations of HAT-P-36 

One transit of HAT-P-36b was observed on April 2013 
through a Gunn-r filter with the Bologna Faint Object 
Spectrograph & Camera (BFOSC) imager mounted on the 
1.52 m Cassini Telescope at the Astronomical Observatory 
of Bologna in Loiano, Italy. The CCD was used unbinned, 
giving a plate scale of 0.58 arcsec pixel -1 , for a total FOV of 
13 arcmin x 12.6 arcmin. Two successive transits were ob¬ 
served on April 2014 by using the Zeiss 1.23 m telescope at 
the German-Spanisli Astronomical Center at Calar Alto, 
Spain. This telescope is equipped with the DLR-MKIII 
camera, which has 4000 x 4000 pixels, a plate scale of 0.32 
arcsec pixel -1 and a FOV of 21.5 arcmin x 21.5 arcmin. 
The observations were performed remotely, a Cousins-/ fil¬ 
ter was adopted and the CCD was used unbinned, too. 
Details of the observations are reported in Tab led As 
in previous uses of the two telescopes (e.g, iMancini et al.l 
l2013 af). the defocussing technique was adopted in all the ob¬ 
servations to improve the quality of the photometric data 
remarkably. Telescopes were also autoguided. 

The data were reduced with a modified version of the 
DEF0T pip eline, written in ID I p. wh ich is exhaustively de¬ 
scribed in ISouthworth et al.l (120141) . In brief, we created 
master calibration frames by median-combining individual 
calibration images and used them to correct the scientific 
images. A reference image was selected in each dataset 
and used to correct pointing variations. The target and 
a suitable set of non-variable comparison stars were iden¬ 
tified in the images and three rings were placed interac¬ 
tively around them; the aperture radii were chosen based 

2 The acronym IDL stands for Interactive Data Language and 
is a trademark of ITT Visual Information Solutions. 


on the lowest scatter achieved when compared with a fit¬ 
ted model. Differential photometry was finally measured 
using the APER routine^. Light curves were created with a 
second-order polynomial fitted to the out-of-transit data. 
The comparison star weights and polynomial coefficients 
were fitted simultaneously in order to minimise the scatter 
outside transit. Final light curves are plotted in Fig. [3] to¬ 
gether with one, obtained by phasing/binning the four light 
curves reported in the discovery paper (IBakos et al.i|2l)12 i. 
shown here just for comparison. 

Very interestingly, the light curves present anomalies 
attributable to the passage of the planetary shadow over 
starspots or starspot complexes on the photosphere of 
the host star, which show up in weak and strong fash¬ 
ion. In particular, the transit observed with the Cassini 
telescope was also obser ved by the Toruri 0.6 m telescope 
dMacieiewski et all 1 20 1 .'ill . which confirms the presence of 
anomalies in the light curves. 

2.3. Photometric follow-up observations of 
WASP-ll/HAT-P-10 

A transit observation of WASP-ll/HAT-P-10b was per¬ 
formed on November 2009 with the Cassini 1.52 m tele¬ 
scope through a Gunn-r filter. Unfortunately, the weather 
conditions were not optimal and the data were badly af¬ 
fected by the poor transparency of the sky. The same tran¬ 
sit monitored by HARPS-N (October 2, 2014) was also ob¬ 
served by the Zeiss 1.23 m telescope through a Cousins- 
/ filter and a more adequate seeing (~ 0.7"). The char- 


3 APER is part of the ASTR0LIB subroutine library distributed 
by NASA. 
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Fig. 4. Light curves of WASP-ll/HAT-P-10b transits compared with the best jktebop fits. The dates, the telescopes 
and the filters used for each transit event are indicated. Residuals from the fits are displayed in t he right panel. St arting 
from the top, the first three light curves are from this work, whereas the last two are from iBakos et all (|2009l ). The 
second and third light curves are related to the same transit event. 


Table 2. Details of the transit observations presented in this work. 


Telescope 

Date of 
first obs 

Start time 
(UT) 

End time 
(UT) 

Nobs 

Texp 

(s) 

T 0 bs 

(8) 

Filter 

Airmass 

Moon 

ilium. 

Aperture 
radii (px) 

Scatter 

(mmag) 

HAT-P-36: 

Cassini 

2013 04 14 

19:58 

00:39 

112 

120 

144 

Gunn r 

1.05 -t 1.000 -> 1.17 

12% 

16,45,65 

1.12 

CA 1.23 m 

2014 04 14 

19:50 

02:08 

176 

80-120 

92-132 

Cousins I 

1.28 1.009 1.20 

98% 

22,50,75 

1.10 

CA 1.23 m 

2014 04 18 

19:54 

00:26 

146 

80-130 

92-142 

Cousins I 

1.22 1.009 1.06 

92% 

22,50,75 

1.30 

WASP-ll/HAT-P-10: 

Cassini 2009 11 14 

23:02 

03:36 

101 

90-150 

114-174 

Gunn r 

1.06 1.02 2.00 

9% 

17,25,45 

2.24 

CA 1.23 m 

2014 10 02 

00:51 

05:23 

125 

120 

132 

Cousins I 

1.12 -)• 1.04 -)• 1.21 

52% 

23,33,60 

0.65 

IAC 80 cm 

2014 10 01 

23:56 

05:39 

179 

90 

115 

Cousins R 

1.008 1.000 2.49 

52% 

28,60,90 

0.92 


Notes. N 0 bs is the number of observations, T exp is the exposure time, T 0 bs is the observational cadence, and ‘Moon ilium.’ is the 
geocentric fractional illumination of the Moon at midnight (UT). The aperture sizes are the radii of the software apertures for the 
star, inner sky and outer sky, respectively. Scatter is the rms scatter of the data versus a fitted model. 


acteristics of the two telescopes and observation modus 
operandi were already described in the previous section. 
Moreover, the identical transit was observed with a third 
telescope, the IAC80cm telescope, located at the Teide 
Observatory on the island of Tenerife (Spain). The opti¬ 
cal channel of this telescope is equipped with CAMELOT, 
a 2048 x 2048 CCD detector with a plate scale of 0.304 arc- 
sec pixel -1 , providing an astronomical FOV of 10.4arcmin 
x 10.4arcmin. Details of these three WASP-ll/HAT-P- 
10 b transit observations are reported in Tabled The data 
sets were reduced in the same way as those for the HAT- 
P-36 case ( Sect. 12.21) and the corresponding light curves 
are shown i n Fig .ffl together with other two taken from 
iBakos et al.l (120091) 7 The simultaneous transit observation 
is highlighted in Fig.jSJ and it allowed an extremely pre¬ 
cise measurement of the mid-point transit time, i.e. To = 
2 456 933.615069 ± 0.000075 BJD(TDB). 


3. Light-curve analysis 


The three HAT-P-36 light curves show possible starspot 
crossing events, which must be analysed using a self- 
consist ent and physically realistic model. As in previous 
cases (|Mancini et al.1 l2013bl l2014bl : iMohler-Fischer et, al.l 
l20l!. we utilise t he PRlShfl and GEMC0 codes 
( Tregloan-Reed et al.l 120131 120151) to undertake this 


task. PRISM performs a modelling of planetary-transit 
light curves with one or more starspots by means of a 
pixellation approach in Cartesian coordinates. GEMC uses 
a Differential Evolution Markov Chain Monte Carlo (DE- 
MCMC) approach to locate the parameters of the PRISM 
model that better fit the data, using a global search. The 
fitted parameters of PRISM are the sum and ratio of the 


4 Planetary Retrospective Integrated Star-spot Model. 

5 Genetic Evolution Markov Chain. 


5 













































































L. Mancini et al.: Physical properties of HAT-P-36b & WASP-ll/HAT-P-10b 


1.005 



0.970 1.. 

933.50 933.55 933.60 933.65 933.70 

BJD(TDB)-2456000 


Fig. 5. Light curve of a transit event of WASP-ll/HAT-P-10b observed simultaneously with the CA 1.23m (red points) 
and the IAC 80 cm (empty orange circles) telescopes. The data of the tw o telescopes are superimposed to highlight the 
difference in the light-curve shape caused by the different filter adopted dKnutson et all 1200711 . The dashed vertical line 
represents the time of transit mid-point that, thanks to the two-site strategy, was very accurately measured. 


fractional radi@ ( r\ + ?’b and k = ^), the orbital period 
and inclination (P and i), the time of transit midpoint 
(To) and the coefficients of the quadratic limb darkening 
law (ua and i>a). 


The prism+GEMC code works in the way that the user 
has to decide how many starspots have to be fitted. We 
fitted for all the spot anomalies that can be seen in the 
transit light curves. Each starspot is then reproduced by 
the longitude and co-latitude of its centre (6 and </>), its 
angular radius (r spo t) and its contrast (p S pot), the last being 
the ratio of the surface brightness of the starspot to that of 
the surroundi ng photosphere. The orbital eccentricity was 
fixed to zero dBakos et al.ll2012lh 


Table 6. Times of transit midpoint of HAT-P-36b and 
their residuals. 


Time of minimum 
BJD(TDB)—2400000 

Cycle 

no. 

O-C 

(JD) 

Reference 

55555.89060 ± 0.00043 

-7 

0.000376 

1 

55597.04088 ± 0.00770 

24 

0.002904 

1 

55601.01882 ± 0.00074 

27 

-0.001197 

1 

55608.98390 ± 0.00030 

33 

-0.000198 

1 

56365.56800 ± 0.00234 

603 

-0.001928 

3 

56397.42170 ± 0.00202 

627 

-0.004655 

3 

56397.42894 ± 0.00031 

627 

0.000822 

2 

56762.44834 ± 0.00018 

902 

-0.000158 

3 

56766.43055 ± 0.00028 

905 

-0.000011 

3 


Notes. References: (1) FLWO 1.2 m ( iBakos et al . 2012 1: (2) 
Cassini 1.52 m (this work); (3) CA 1.23m (this work). 


The derived parameters of the planetary system are re¬ 
ported in Tabled while those of the starspots in Tabled 
The light curves and their best-fitting models are shown in 
Fig. [3] For each transit, a representation of the starspots on 
the stellar disc is drawn in Fig. [ 5 ] 


The light curves for the WAS P-ll/HAT-P -10 transits 
were modelled by jkteboiB (see ISouthworthi [2008 . 120131 
and references therein), as this code is much faster than the 
previous one and no clear starspot anomalies are visible. 
The parameters used for JKTEBOP were the same as for 
PRISM and the orbital eccentricity was also fixed to zero 
(Bak os et al. 20091) . Ligh t curves from the second discovery 
paper dBakos et al.ll2009 f were also reanalysed. The results 
of the fits are summarised in Tableland displayed in Fig.0] 
Its low quality meant that the Cassini light curve (top curve 
in Fig. H was excluded from the analysis. 


3.1. Orbital period determination 

We used the new photometric data to refine the transit 
ephemeris of HAT-P-36b. The transit times and uncertain¬ 
ties were obtained using PRISM+GEMC as previously ex- 
plaine d. To these timings , we added four from the discovery 
paper jBakos et, al.ll2012f) and two from iMacieiewski et, al.l 
(l2013ln . The nine timings were placed on the BJD(TDB) 
time system (Tabled. The resulting measurements of tran¬ 
sit midpoints were fitted with a straight line to obtain a 
final orbital ephemeris 

P = 1.32734684 ± 0.00000050 d 

6 The fractional radii are defined as ta = — and n, = —, 
where A* and R p are the true radii of the star and planet, and 
a is the orbital semimajor axis. 

' jktebop is written in FORTRAN77 and is available at: 
http://www.astro.keele.ac.uk/jkt/codes/jktebop.html 
s We added 1/25 days to the published values to correct an 
error caused by a misunderstanding in reading the time stamps 
of the fits files. 
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Table 3. Photometric properties of the HAT-P-36 system derived by fitting the light curves with GEMC. 


Source 

Filter 

va + Tb 

k 

i° 

UA 

VA 

Cassini 1.52 m 
CA 1.23 m #1 
CA 1.23 m #2 

Gunn r 
Cousins I 
Cousins I 

0.2278 ±0.0037 
0.2298 ±0.0039 
0.2285 ±0.0039 

0.1243 ±0.0025 
0.1290 ±0.0023 
0.1298 ±0.0012 

85.95 ±0.34 
85.76 ± 0.36 
85.86 ± 0.43 

0.24 ±0.13 
0.25 ±0.10 
0.32 ±0.11 

0.40 ±0.21 
0.29 ±0.19 
0.15 ±0.14 


Table 4. Starspot parameters derived from the GEMC fitting of the HAT-P-36 transit light curves presented in this work. 


Telescope 

starspot 

0(°) a 

M°) b 

r spo t(°) c 

d 

P spot 

Temperature (K) e 

Cassini 1.52 m 

#1 

#2 

-36.59 ± 1.85 
38.29 ± 2.88 

77.41 ± 5.82 
82.68 ± 5.27 

11.42 ±3.45 
12.54 ±5.57 

0.65 ±0.07 
0.62 ±0.23 

5070 ± 129 

5016 ±421 

CA 1.23 m #1 

#1 

21.06 ±3.20 

86.52 ±3.96 

17.67 ±4.81 

0.86/°;“ 

5380 ± 305 

CA 1.23 m #2 

#1 

#2 

-24.99 ±4.91 
46.00 ± 14.13 

87.08 ±3.17 
88.27 ±2.38 

21.38 ±6.82 
18.66 ±3.52 

0.92/°'°! 
0.78 ±0.19 

5485 ± 244 

5253 ± 353 


Notes. ^The longitude of the centre of the spot is defined to be 0° at the centre of the stellar disc and can vary from —90° 
to 90°. ^The co-latitude of the centre of the spot is defined to be 0° at the north pole and 180° at the south pole. ^Angular 
radius of the starspot (note that an angular radius of 90° covers half of stellar surface). ^Spot contrast; note that 1.0 equals the 
brightness of the surrounding photosphere. *- e) The temperature of the starspots are obtained by considering the photosphere and 
the starspots as black bodies. 


Table 5. Photometric properties of the WASP-ll/HAT-P-10 system derived by fitting the light curves with jktebop. 


Source 

Filter 

ta ± r h 

k 

i° 

UA 

VA 

Cassini 1.52 m 

Gunn r 

0.10137 ±0.01380 

0.13639 ±0.00868 

88.00 ± 1.97 

0.56 ±0.23 

0.21 ±0.06 

CA 1.23 m 

Cousins / 

0.09230 ± 0.00197 

0.13234 ±0.00078 

89.09 ±0.68 

0.24 ±0.04 

0.25 ±0.07 

IAC 80 cm 

Cousins R 

0.09695 ± 0.00358 

0.13657 ±0.00200 

88.29 ± 0.54 

0.40 ±0.08 

0.21 ±0.07 

FLWO 1.2 m 

Sloan z 

0.09277 ± 0.00350 

0.12818 ±0.00181 

89.95 ±0.92 

0.41 ±0.06 

0.27 ±0.07 

FLWO 1.2 m 

Sloan i 

0.09160 ± 0.00242 

0.12856 ±0.00156 

89.94 ± 0.76 

0.49 ±0.05 

0.25 ±0.07 


and the mid-transit time at cycle zero 

T 0 = 2455565.18165 ± 0.00037 BJD(TDB), 

with reduced xt = 2-5, which is significantly greater than 
unity. A plot of the residuals around the fit (see Fig. EH) 
does not indicate any clear systematic deviation from the 
predicted transit times. 

We introduced other 21 timings in the analysis mea¬ 
sured based on transit light curves observed by amateur 
astronomers and available on the ETE0 website. These 21 
light curves were selected considering if they had complete 
coverage of the transit and a Data Quality index < 3. 
Repeating the analysis with a larger sample, we obtained 
4.1 as the reduced xt of the fit. A high value of xl means 
that the assumption of a constant orbital period does not 
agree with the timing measurements. This could indicate 
either transit timing variations (TTVs) or that the mea¬ 
surement errors are underestimated. The latter possibility 
can occur very easily, but be difficult to rule out, so the de¬ 
tection of TTVs requires more than just an excess xt- The 
Lomb-Scargle periodogram generated for timing residuals 
shows no significant signal, making any periodic variation 
unlikely, so we do not claim the existence of TTVs in this 
system. The residuals of the ETD timings are also shown 
in Fig. [6] for completeness. 

New mid-transit times (three transits; this work) and 
those available in the literature (9 transits) allowed us to 

9 The Exoplanet Transit Database (ETD) website can be 
found at http://var2.astro.cz/ETD 


Table 7. Times of transit midpoint of WASP-ll/HAT-P- 
10 b and their residuals. 


Time of minimum 
BJD(TDB)—2400000 

Cycle 

no. 

O-C Reference 
(JD) 

54729.90657 ± 0.00023 

0 

-0.000589 

1 

54770.85432 ±0.00011 

11 

-0.000115 

1 

54793.19692 ± 0.00076 

17 

0.002071 

2 

55143.10895 ± 0.00046 

111 

0.001013 

2 

55161.71529 ±0.00021 

116 

0.000490 

3 

55507.90419 ± 0.00042 

209 

-0.001219 

3 

55842.92921 ± 0.00021 

299 

0.000631 

3 

55842.92952 ± 0.00044 

299 

0.000941 

3 

55865.26734 ± 0.00166 

305 

-0.001483 

2 

55150.54662 ± 0.00104 

113 

-0.000741 

4 

56933.61504 ± 0.00008 

592 

-0.000086 

5 

56933.61522 ± 0.00018 

592 

0.000099 

6 


Notes. References: (1) FLWO 1.2 m llBakos et, ahl l200fi)i : (2) 
iWang et all d2014f >: (3) iSada et ahl (120121) : 141 Cassini 1.52 m 
(this work); (5) CA 1.23m (this work); (6) IAC 80cm (this 
work). 


refine transit ephemeris for WASP-ll/HAT-P-10 b tooQn 
particular, the two light curves from iBakos et alJ (12009} ) 
were re-fitted with jktebop (Fig. [I]). All the timings are 
summarised in TablejT] As a result of a linear fit, in which 
timing uncertainties were taken as weights, we derived 

P = 3.72247967 ± 0.00000045 d 
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Fig^. Resi duals for the timings of HA T-P-3 6 b at mi d-tran sit versus a linear ephemeris. The four timings reported by 
iBakos et al.l (120121 ) and the two by iMacieiewski et ahl (120131) are plotted using empty circles, while those marked with 
gray circles are from ETD. The timings measured by this work are shown with black circles. The dashed lines mark 2-cr 
uncertainties of the ephemeris (i.e. at the 95% confidence level), which was estimated without considering the data from 
ETD. 


and the mid-transit time at cycle zero 

T 0 = 2454729.90716 ± 0.00020 BJD T db , 

with xi = 5-0, which is again far from unity. This could be 
explained again by underestimated timing uncertainties or 
hypothesizing a variation in transit times caused by unseen 
planetary companion or stellar activity. However, also for 
this case the Lomb-Scargle periodogram shows no signifi¬ 
cant signal, not supporting the second hypothesis. The plot 
of the residuals around the fit is shown in Fig. [7] 

As in the previous case, we selected 34 transits from 
ETD, using the same criteria explained above. Most of the 
corresponding timings are very scattered around the pre¬ 
dicted transit mid-times, but have tiny error bars. Including 
them in the fit returns a higher value for the suggesting 
that their uncertainties are very underestimated, and the 
data reported on ETD should be used cautiously. They are 
shown in Fig. [7] 

We would like to stress however that, since in both the 
cases the timing measurements are in groups separated by 
several hundred of days, we are not sensitive to all the 
periodicities. Only by performing systematic observations 
of many subsequent transits it is possible to rule out the 
presence of additional bodies in the two TEP systems with 
higher confidence. 


3.2. HAT-P-36 starspots 

As described in Sect.[3j the anomalies in the three HAT-P- 
36 light curves were modelled as starspots, whose param¬ 
eters were fitted together with those of the transits (see 
Tables [3] and S]) . The stellar disc, the positions of starspots 
and the transit chords are displayed in Fig. [51 based on 
the results of the modelling. Considerin g bo th the photo¬ 
sphere and the starspots a s black bodies ( RabuseLaL 20091: 
ISanchis-Oieda et ahl 20111 iMohler-Fischer et al. 12013 ) and 


oancnis-wiecia ei ai. zun; iviomcr-r iscucr cl ai. zuio, 

using Eq. l' of IsITral (f2003il and Tbr = 5620 ± 40 

Sect. 14.II) . we estimated the temperature of the stars] 


at different bands and reported them in the last column of 
Table [IJ The values of the starspot temperature estimated 
in the transit on April 2013 are in good agreement with 
each other within the experimental uncertainties. The same 
is true for those observed in the two transits on April 2014, 
even if they point to starspots with temperature higher 
than, but still compatible, with those of the previous year. 
The starspot temperatures measured from the Loiano and 
CAHA observations are consistent with what has been mea¬ 
sured for other main-sequence stars in transit observations, 
as we can see from Fig. [9] where we report the starspot 
temperature contrast versus the temperature of the photo¬ 
sphere of the corresponding star for data taken from the lit¬ 
erature. The spectral class of the stars is also reported and 
allows us to figure out that the temperature difference be¬ 
tween photosphere and starspots is n ot str ongly de pende nt 
on spectral type, as already noted by IStrassmeierl (120091 1. 

The observations of multiple planetary transits across 
the same starspot or starspot complex could provide an¬ 
other type of precious information dSanchis-Oieda et al.l 
1201 If) . Indeed, thanks to the good alignment between the 
stellar spin axis and the perpendicular to the planet’s or¬ 
bital plane, one can measure the shift in position of the 
starspot between the transit events and constrain the align¬ 
ment between the orbital axis of the planet and the spin axis 
of the star with higher p recisio n than from the measurement 
of the RM effect (e.g.. lTregloan-Reed et al.ll2013f l. 

In our case, we have two consecutive transits observed 
in 2014 April 14 and 18 and we might wonder if the planet 
has crossed the same starspo t complex in t hose tr ansit 
events. According to Eq. 1 of iMancini et al.l (|2014bf ). the 
same starspot can be observed after consecutive transits or 
after some orbital cycles, presuming that in the latter case 
the star performs one or more rotations around its axis. 
Since the projected obliquity measured by the RM effect 
(see Sect.[4]) indicates spin-orbit alignment, we expect to 
find similar values for the starspot parameters in our fits. 
Examining Tabled we note that they seem to agree on co- 
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Fig. 7. Residuals for the timings of WASP-ll/HAT-P-l Ob at mid-transit versus a linear ephemeris. The filled black dots 
are from this work, open circles from the literature, (i.e. iBakos et al.ll2009t ISada et aI1l2012t IWang et al.ll2014T) and gray 
dots from ETD. The dashed lines mark 2-a uncertainties of the ephemeris, which was estimated without considering the 
data from ETD. 




Fig. 8. Representation of the stellar disc, starspot positions and transit chord for the three HAT-P-36 datasets containing 
spot anomalies: (a) refers to the transit observed with the Cassini 1.52 m telescope on 2013 April 15, while (b) and (c) 
to those observed with the CA 1.23 m telescope on 2014 April 14 and 18, respectively. The gray scale of each starspot is 
related to its contrast (compare with values in Tablc|4|). 


latitude, size and contrast within their 1-a errors. Owing 
to its size, the starspot on April 14 should still be seen 
four days later (if it is on the visible side of the star). That 
means that the starspot should have travelled ~ 315° in 
four days, giving a rotation period of ~ 4.57 days, which 
is too fast compared with that at the stellar equator esti¬ 
mated from the sky-projected rotation rate and the stellar 
radius, that is 

Piot ~ sini* = (16.9 ± 4.3 d) sini*, (1) 

v sin i * 

where i* is the inclination of the stellar rotation axis with 
respect to the line of sight. If we consider reversing the 
rotational direction so that the starspot is travelling right 
to left, then the starspot had to travel for ~ 45° in four 
days; but this implies that the stellar rotation period at 
co-latitude of 87° would be now too large (~ 32 days) and 
that HAT-P-36 b has a retrograde orbit. Since the latter 


hypothesis is excluded by the geometry of the RM effect 
that we observed (see Sect. 14.3D and since four days are not 
sufficient for the starspot to rotate around the back of the 
star and then appear on the left hemisphere, we conclude 
that the starspot observed on 2014 April 14 is different from 
those on April 18. 

3.3. Frequency analysis of the time-series light curves 

In addition to the new photometric measurements of 
the transits presented in this work, dense time-series 
light curves are available in the WASP and HAT 
archives/databases for both the stars. These data are po¬ 
tentially very useful for detecting the signatures of the ro¬ 
tational period of each of the two stars. This is particularly 
true for the case of HAT-P-36, since the starspots detected 
in the photometric light curves (see previous section) indi¬ 
cate that this star is active. A further confirmation of this 
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Fig. 9. Starspot temperature contrast with respect to the photospheric temperature in several dwarf stars. The name 
and spectral type of the st ar are also rep orted. The values for TrE S-1, CoRoT-2, HD 189733 , WASP-4, HATS-2, Kepler - 
63, Qatar-2 are taken from iRabus e t al. (20 09lf. Silva -Valio et al.1 (1 201 f ). ISimr et, al.l (1201 ill . ISanchis-Oieda et all (|201lh . 
iMohler-Fischer et al.1 (1201 311 .ISanchis-Oieda et al.l (te oTshlMancini et a" (l2014blf . respectively. The t wo values for WASP - 
19 are from Mancini et al.1 (l2013blf and lHuitson et al. ( 201311 . The value of the Sun was taken from iBerdvueinal (I2005H . 
The error bars have been suppressed for clarity. Note that some stars appear twice or more. The black dots refer to 
the values estimated for the case of HAT-P-36 (this work). In particular, the lower is the value estimated by taking the 
weighted mean from the starpots detected with the Cassini telescope in April 2013, while the higher value was calculated 
from starpots detected with the CA 1.23 m telescope in April 2014. 


activity is also given by the spectral analysis, as seen in 
Sect. 14.21 

We retrieved the photometric time series for HAT-P- 
36 (SDSS r filter) and WASP-ll/HAT-P-10 (/ Bessel fil¬ 
ter) from the HAT public archivcE! more s pecific ally , we 
used the magnitudes data sets tagged as TF1 (iKovacs et alJ 
2(lOA . The WASP data available for WASP-ll/HAT-P-10 
were downloaded from the NASA Exoplanet ArchivcEl. The 
in-transit data points were removed from the light curve of 
each dataset. 

The frequency analysis of the WASP and HAT time- 
series data related to WASP-ll/HAT-P-10 did not detect 
any significant periodicity. Only the HAT data show a non¬ 
significant peak at / = 0.034 d -1 . It is close to the synodic 
month, and the folded light curve shows a gap in the phase 
coverage corresponding to the Full Moon epochs. Therefore, 
it has been recognized as a spurious peak due to a combi¬ 
nation of instrumental effects and spectral window. 

On the other hand, the analysis of the HAT time-series 
light curve of HAT-P-36 (see Fig. [10]) provided more inter- 


10 http://hatnet.org 

11 http://exoplanetarchive.ipac.caltech.edu/ 


esting results. After removing the transit data, the time- 
series light curve is composed of 9891 measurements. Both 
the iterative sine-wave fitting method (ISWF; lYanicekl 
1197111 and the Lomb-Scargle periodogram supplied a lower 
peak at f = 0.066 d -1 and a higher one at 2/ = 0.13 d -1 . 
Figured]] shows the power spectrum obtained with the 
ISWF method. The amplitudes of the 2/ and / components 
are 1.40±0.15 and 0.88±0.15 mmag, respectively. Since the 
level of the noise is 0.2 mmag, bot h are greater than th e 
S/N — 4.0 significance threshold (iKuschnig et al.l H 9071) . 
Moreover, the harmonics 3/ and 4/ are both clearly visible, 
the latter being more relevant than the former. HAT-P-36 
was measured several times per night and often over night. 
Hence we could calculate 97 night averages (red points in 
Fig. El. The frequency analysis of these averages supplied 
the same results as the single measurements. 

When the even harmonics have amplitudes larger than 
those of the odd ones, the resulting light curve is shaped as 
a double wave. Since we know that the star is seen equator- 
on (see Sect. S3D, such a light curve can easily be due to 
two groups of starspots that are alternatively visible to the 
observer. This is an independent confirmation of the spot 
model suggested in the previous subsection. We determined 
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Phase 


Fig. 10. Top panel: Single measurements (grey dots) and 
average values on single nights (red circles with errorbars) 
folded with P = 15.6 d, the best period obtained fitting 
the data with / and 2/ fixed. Bottom panel: Single mea¬ 
surements (grey dots) and average values (red circles with 
errorbars) folded with P = 15.2 d, the best period obtained 
leaving the frequencies free to vary. 


the rotational period in two ways. First, we searched for the 
best fit by fixing the frequency / and its harmonic 2/. This 
procedure yields P ro t = 15.2 ± 0.3 d. Then, we supposed 
that the differential rotation and the spread in longitude 
did not justify this condition, so we left the two frequency 
values free to vary, obtaining P rot = 15.6 ± 0.1 d. The pe¬ 
riods are the same within the error bars and in excellent 
agreement with the period inferred from the stellar param¬ 
eters, see Eq. ©• We merged the two determinations into 
Prot = 15-3 ± 0.4 d to cover both values and respective er¬ 
rorbars. The HAT-P-36 light curve is shown in Fig. flT)l 
phase-folded with P rot = 15.6 (top panel) and P ro t = 15.2 
(bottom panel) days. 



Fig.11. Power spectrum of the original data of HAT-P-36, 
after removing the measurements during transit events. 

4. HARPS-N spectra analysis 

4.1. Spectroscopic determination of stellar atmospheric 
parameters 

The effective temperature (P e ff), surface gravity (log< 7 *), 
iron abundance ([Fe/H]) and microturbulence velocity 
(i’mic) of the two host stars we re derived by using the spec¬ 
tral analysis package MOOG (ISnedenlll973l version 2013), 
and the equiv alen t widths (EW s) of iron lines, as described 
in detail in lBiazzo et al . 1 (j2012f h The EWs were measured on 
the mean spectrum obtained by averaging all the HARPS-N 
spectra available for each of the two stars. Here, T e g was de¬ 
termined by imposing the condition that the Fe I abundance 
does not depend on the excitation potential of the lines, the 
microturbulence velocity, by requiring that the Fel abun¬ 
dance is independent on the line EWs, and logg* by the 
Fel/Fen ionization equilibrium condition. The results are 
reported in Tables[5]and[9l for HAT-P-36 and WASP-11, re¬ 
spectively. The projected rotational velocity was estimated 
by means of spectral synthesis analysis with MOOG, yield¬ 
ing usini* values of 3.0 ± 0.5 kms" 1 and 0.9 ± 0.5 kms" 1 
for HAT-P-36 and WASP-11, respectively, which are con¬ 
sistent, within the errors, with those obtained from the RM 
effect modelling (see Sect. 14.311 . 

4.2. Stellar activity index 

The median values of the stellar activity index estimated 
from the HARPS-N spectra turned out to be logi?(j K = 
—4.636 ± 0.066 dex for HAT-P-36 and -4.848 ± 0.029 dex 
for WASP-ll/HAT-P-10, indicati ng a moderate and low 
activity, respectively le.g. iNoves et al.lfl984lh These values 
agree with what we have found in Sect. 13.31 by analysing 
the big photometric data sets from the HATNet amd 
WASP surveys. The above values were obtained by adopt¬ 
ing B — V = 0.719 and 0.989, respectively, based on the stel¬ 
lar effective temperat u re and the color- T R g co n versio n table 
by ICasagrande et al.l (|2010H . iKnutson et all (|2010fl found 
logi?Jj K = —4.823dex for WASP-ll/HAT-P-10 (estimated 
considering B — V = 1.01), which is in good agreement with 
our value. 

Thanks to the activity index, we can also assess the 
expected rotation period from the level of the stellar activ¬ 
ity, obtain ing for HAT-P-36 P rn t, = 18.3 ± 3.1 d and 18.7 ± 
2.9 d using lNoves et alJ (jl984li and lMamaiek fc Hillenbrandl 
(l2008i l calibration scales, respectively, in quite good agree¬ 
ment, within the observational uncertainties, with what 
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was estimated in Sect. rm For WASP-ll/HAT-P-10, we 
obtained 40.8 ± 1.8 d and 40.9 ± 2.2 d, respectively. 

Finally, we can get a clue on stellar age by 
app lying the acti vity-age calibration proposed by 
iMamaiek fc Hillenbrandl (|2008f) , obtaining 2 Gyr and 
5Gyr for HAT-P-36 and WASP-ll/HAT-P-10, respec¬ 
tively. However, we notice that the ages derived from 
stellar rotation and activity could be altered in the TEP 
systems because of star-planet tidal interactions, and any 
discrepancy with the values derived from the isochrones 
(see Sect. [5]) may be due to this. 

4.3. Determination of the spin-orbit alignment 

The analysis of the photometric data allowed us 
a direct measurement, without reliance on theoreti- 
cal stellar models, of the stellar mean densit y, p* 
(|Seager fc Mallen-Ornelasl l200.il ISozzetti et al.112007H and, 
in combination with the spect roscopic orbital solution , the 
planetary surface gravity, g p (|Southworth et all 120071) . for 
each system. Exploiting the measured values of p*, T e g and 
[Fe/H], we made use of the Yonsei-Yale evolutionary tracks 
ijDemarciue et al.lfeOOdH to determine stellar characteristics, 
including M* and f?*. The HARPS-N RV data sets were 
then fitted using a model that accounts both for the RV 
orbital trend and the RM anomaly. 

We us ed the RM model t hat was elaborated a nd de - 
scribed in ICovino et al.1 (1201311 and lEsposito et al.1 (I2014H , 
and we used the same least-square minimization algorithm 
to adjust it to the data. We set as free the barycentric radial 
velocity, 7 , the mass of the planet, M p , the projected stellar 
rotational velocity, v sini*, the projected spin-orbit angle, 
A and the linear limb-darkening coefficient u. All the other 
relevant parameters were kept fixed to the values deter¬ 
mined from the spectroscopic and light curves analysis. The 
best-fitting RV models are illustrated in Fig.|T]and[2l super¬ 
imposed on the datasets. In particular, the sky-projected 
spin-orbit misalignment angle was A = —14 ± 18 deg for 
HAT-P-36 and A = 7 ± 5 deg for WASP-ll/HAT-P-10, in¬ 
dicating an alignment of the stellar spin with the orbit of 
the planet for both the systems. The final error bars were 
computed by a bootstrapping approach. The greater uncer¬ 
tainty in HAT-P-36 is of course because the transit was not 
completely observed. The reduced iQ maps in the A— v sin i* 
plane do not point out any clear correlation between A and 
v sini*. 

The precise knowledge of the rotational period of a star, 
allows estimating the stellar spin inclination angle, «*, and 
thus the true spin-orbit alignment angle, ip. This is exactly 
the case for HAT-P-36, for which the analysis of the long 
time-series data set recorded by HATNet survey allowed us 
to determine the stellar rotational period. By using P rot = 
15.3 ± 0.4d (Sect. 13.31) . we estimated that u = 65° ± 34°. 
Then, we used Eq. (7) in lWinn et al.l (l2007lh 

cos ip = cos i* cos i + sin i* sin i cos A, (2) 

to derive the true misalignment angle, obtaining ip = 25 I 25 . 
The major source of error in the determination of both ?’* 
and ip, estimated through the propagation of uncertainties, 
comes from the large relative error in the measurement of 
csin**; the observation of a full transit would have reduced 
the errors by a factor of 2 or more. 


5. Physical parameters of the two systems 


For calculating the full physical properties of each pla netary 
syste m, we used the HSTEP methodology (see lSouthworthl 
120121 . and references therein). The values of the photomet¬ 
ric parameters rx + n,, k and i were first combined into 
weighted means. The orbital eccentricity was fixed to zero. 
We added in the measured P from the light curve analysis, 
T e g and [Fe/H] from the spectroscopic analysis, and the ve¬ 
locity semi-amplitude of the star, K a, measured fr o m the 
RVs . For Ka we used the values from iBakos et al.l ((2012) 
and lBakos et al.l (2009). These come from spectra with rea¬ 
sonable coverage of the orbits of the host stars, so are more 
precise than our own values which come from data obtained 
only during or close to transit. We estimated a starting 
value for the velocity semi-amplitude of the planet, K p , and 
used this to determine a provisional set of physical prop ¬ 
erties of the system using standard formulae (e.g.. iHilditchl 

(MU). 


We then interpolated within a set of tabulated predic¬ 
tions from a theoretical stellar model to find the expected 
stellar radius and T e ff for our provisional mass and the ob¬ 
served [Fe/H]. The value of K p was then iteratively refined 
to maximise the agreement between the observed and pre¬ 
dicted T e ff, and the provisional (Rx + Rb)/a and measured 
rx + H> values. This Kb was then used alongside the other 
quantities given above to determine the physical properties 
of the system. This process was performed over all possible 
ages for the stars, from the zero-age to the terminal-age 
main sequence, and the overall best-fitting age and result¬ 
ing physical parameters found. 

The above process was performed usi ng each of five dif- 
ferent sets of theoretica l models: Glaret dG laret 2 004(1 . Y 2 
(iDemarq ue et al.| 2004(1 . BaSTI dPietrinferni et ah 2004 b 
VRS S ( VandenBerg et ail 1200611 and DSEP ( Potter et al l 
12008(1 . The final result of this analysis was five sets of prop¬ 
erties for each system, one from each of the five sources 
of theoretical models (see Tables lATl and IA.2I) . For the fi¬ 
nal value of each output parameter, we calculated the un¬ 
weighted mean of the five estimates from using the different 
sets of model predictions. Statistical errorbars were prop¬ 
agated from the errorbars in the values of all input pa¬ 
rameters. Systematic errorbars were assigned based on the 
inter-agreement between the results from the five different 
stellar models. 

Tables [ 8 ] and [9] give the physical properties we found for 
HAT-P-36 and WASP-11, respectively. Our results are in 
good agreement, but they are more precise than previous 
determinations. 


6. Discussion 

As introduced in Sect.[l] orbital obliquity could be an im¬ 
portant parameter to determine in order to trace the phys¬ 
ical processes that happened during migration phase of 
current hot Jupiters. Several empirical trends have been 
presented to corroborate this hypothesis. TEP systems, in 
which the parent star is relatively cool (T e ff 6100 K), 
should be much more aligned than those with hotter stars, 
because the realignment of cool star’s convective envelopes 
can occur on timescales shorter than planet’s orbital deca y 
((Winn et, al.l l201fl( lAlbrecht et al.1 l2012bl : |D awso Til 12014ft . 
The larger spin-orbit misalignment expected for the hot 
stars (F-type stars with T e s > 6100 — 6200 K) could be 
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Table 8. Physical parameters of the planetary system HAT-P-36 derived in this work. 


Parameter 

Nomen. 

Unit 

This Work 

Bakos et al. (2012j 

Stellar parameters 

Spectral class . 

Effective temperature . 

Metal abundance . 

Projected rotational velocity 

Rotational period . 

Linear LD coefficient . 

Mass . 

Radius . 

Mean density. 

Logarithmic surface gravity . 
Age . 

Teff 
[Fe/H] 
v sini* 
Plot 
u 

M* 

R* 

P* 

log 3* 

K 

kms" 1 
days 

M 0 

R & 

PQ 

cgs 

Gyr 

G5V 

5620 ± 40 
+0.25 ±0.09 

3.12 ±0.75 

15.3 ±0.4 

0.81 ±0.17 

1.030 ±0.029 ±0.030 

1.041 ±0.013 ±0.010 

0.913 ±0.027 

4.416 ±0.010 ±0.004 

A c+2.4+3.1 
-2.8 

5560 ± 100 
±0.26 ±0.10 

3.58 ± 0.5 

1.022 ± 0.049 

1.096 ± 0.056 

4.37 ±0.04 

6-6 ±f;| 

Planetary parameters 

Mass . 

M p 

Mj up 

1.852 ±0.088 ±0.036 

1.832 ± 0.099 

Radius . 

R P 

Rjup 

1.304 ±0.021 ±0.013 

1.264 ±0.071 

Mean density. 

Pp 

PJup 

0.737 ± 0.095 ± 0.007 

0.84 ±0.14 

Surface gravity . 

9p 

ms -2 

27.0 ±1.4 

28 ±3 

Equilibrium temperature ... 
Safronov number . 

T eq 

0 

K 

1788 ± 15 

0.0658 ± 0.0030 ± 0.0006 

1823 ± 55 

0.067 ± 0.005 

Orbital parameters 

Time of mid-transit . 

To 

BJD(TDB) 

2 455 565.18167 ±0.00036 

2 455 565.18144 ± 0.00020 

Period . 

Torb 

days 

1.32734683 ± 0.00000048 

1.327347 ± 0.000003 

Semi-major axis . 

a 

au 

0.02388 ± 0.00022 ± 0.00023 

0.0238 ± 0.0004 

Inclination . 

i 

degree 

85.86 ±0.21 

86.0 ± 1.3 

RV-curve semi-amplitude ... 

K a 

ms -1 

316 ± 39 “ 

334.7 ± 14.5 

B ary centric RV. 

7 

km s —1 

-16.327 ±0.006 

-16.29 ±0.10 

Projected spin-orbit angle .. 

A 

degree 

-14 ±18 


True spin-orbit angle . 

ip 

degree 

oc+38 
zo — 25 



Notes. Where two errorbars are given, the first refers to the statistical uncertainties and the second to the systemat ic errors. 
a This value of Ka was determined from our RV data, which were taken during transit time only. The value reported bv lBakos et all 
( l2009h . which is based on RV data with a much better coverage of the orbital phase, was therefore preferred for the determination 
of the other physical parameters of the system (see Sect(5J. 


also explained by considering they have a smaller convec¬ 
tive zone than that of the cool stars. This implies that the 
tidal dissipati on is much less efficient for the hot stars than 
the cold ones (|Laill2012[ IValsecchi fc Rasioll2014lh 

It was also noted that TEP systems in which the hot 
Jupiter is massive (M p > 3 Mj up ) tend to have lower spin- 
orbit angles, the parent star bein g much more affected by 
planet’s tidal influence (iHebrard et al.ll201lh . However, sev¬ 
eral theoretical s tudies were not able to explain the above 
correlations (e. g., Rogers fc Linll2013l ; lXue et al.ll2014f) and, 


as stressed bv fesposito et all (|2014fl . there are planetary 


systems composed by cool stars that are highly misaligned. 

With this work, we contribute to enlarge the sam¬ 
ple, adding the measurement of the sky-projected orbital 
obliquity for other two TEP systems. They are HAT-P- 
36 and WASP-11/HAT-P10, both composed by stars with 
T e s < 6100 K, whose spin results to be well aligned with the 
planetary-orbit axis. They are reported in Fig. [12] (black cir¬ 
cles), together with the other 82 known TEP svstemJ^L as 
a function of stellar effective temperature, projected rota¬ 
tional velocity, eccentricity, orbital perio d, planetary mass 
and stellar metallicity. Following iDawso , we di¬ 

vided the data in two groups, based on the parent stars 
temperature; they are shown with orange circles if T e s > 
6100 K and with empty circles for the opposite case. An in¬ 
spection of these plots does not allow to highlight any clear 


12 The data were taken from TEPcat llSouthworthll20l'ij 'l. The 
brown dwarfs KELT-1 and WASP-30 were excluded from the 
sample. 


trend, but, on the contrary, the values of A are quite ran¬ 
domly distributed in the various parameter spaces, suggest¬ 
ing that the observed diversity of stellar obliquities may be 
a consequence of more complicated interactions with outer 
planets. 

We also plotted in Fig. [13] the projected orbital obliq¬ 
uity o f the cool-star systems, i.e., following I Anderson et al.l 
(1201511 . those with T eff < 6150 K, as a function of orbital 
distance in units of stellar radii, a/1?.*. However, we do not 
see an absolute confinement of A at any particular orbital- 
separation range. 

According to TEPCat, there are now 12 TEP systems 
for which we have the measurement of the true obliquity, i.e. 
the angle, ip, between the axes of the stellar rotation and the 
planetary orbit. These values are plotted against the stellar 
temperature in Fig. [14] that, again, does not highlight any 
correspondence between the two quantities. 

Based on a recent gyrochronology work by 
iMeibom et al.l (l2015f) . the rotation period of HAT-P- 
36, estimated in Sect. Efl should correspond to an age of 
about 1.8 Gyr, which is close to the lower limit reported 
in Table(5] as derived from our isochrone-fitting analysis. 
Since the planet is very close to its host, we may speculate 
about tidal interactions affecting the evolution of the orbit 
and of the stellar spin. Assuming a modified tidal quality 
factor Q' = 10 6 , the remaining lifetime of the system 
- i.e., the time left bef ore the planet plunges into the 
star - is only 35 Myr dMetzger et alJl2012h . Such a rapid 
orbital decay should be accompanied by a negative transit 
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Table 9. Physical parameters of the planetary system WASP-ll/HAT-P-10 derived in this work. 


Parameter 

Nomen. / Unit 

This Work 

West et al. (2009) 

Bakos et al. (2009) 

Wane et al. (2014') 

Stellar parameters 

Spectral Class . 


K3V 

K 

K 


Effective temperature. 

T eff (K) 

4900 ± 65 

4800 ± 100 

4980 ± 60 


Metal abundance . 

[Fe/H] 

+0.12 ± 0.09 

+0.0 ±0.2 

+0.13 + 0.08 


Proj. rotational velocity .. 

usini* (kms x ) 

1.04 + 0.15 

< 6.0 

0.5 ± 0.2 


Linear LD coefficient . 

u 

0.78 + 0.16 




Mass . 

(M 0 ) 

0.806 + 0.038 + 0.013 

0 77+°- 10 
u - 1 ‘ -0.08 

0.83 ± 0.03 

0.862 + 0.014 

Radius . 

-R* (Rq) 

0.772 + 0.014 + 0.004 

0 74+ 0 04 

u -' -0.03 

0.79 ± 0.02 

0.784/!;!!! 

Mean density . 

P* ( PQ ) 

1.748 ± 0.074 



1.789/!;°!! 

Logarithmic surface gravity 

logs* (cgs) 

4.569 + 0.018 + 0.002 

4.45 + 0.20 

4.56 ± 0.02 


Age . 

(Gyr) 

ij £+5.8 +1.6 

1 - u —3.0 -1.8 


7.9 ±3.8 


Planetary parameters 

Mass . 

M p (Mj up ) 

0.492 ± 0.023 ± 0.005 

0.53 + 0.07 

0.487 + 0.018 

0.526 + 0.019 

Radius . 

Rp (Rjup) 

0.990 ± 0.022 ± 0.005 

o-9iio:os 

1 OO^ -1-11,032 
i.uuo_o 027 

f) 000+^.0^9 

u.yyy_Q oi8 

Mean density . 

Pp (pjup) 

0.475 ± 0.026 ± 0.002 

0.69±°;°I 

0.594 ± 0.052 

0 526”*" 0 035 

u.ozu_q Q46 

Surface gravity . 

g p (ms -2 ) 

12.45 + 0.50 

14.45/j;l! 

12.02 + 0.83 


Equilibrium temperature . 

Teq (K) 

992 ± 14 

960 ± 70 

1020 + 17 

1006.5/1°! 

Safronov number . 

e 

0.0539 ± 0.0019 ± 0.0003 


0.053 + 0.002 


Orbital parameters 






Time of mid-transit . 

T 0 (BJD tdb ) 

2 454 729.90716 
±0.00020 

2 454 473.05588 
±0.00020 

2 454 759.68683 
±0.00016 

2 454 808.07904 
±0.00012 

Period. 

Rorb (days) 

3.72247967 

3.722465 

3.7224747 

3.72247669 

±0.00000045 

±0.000007 

±0.0000065 

±0.00000181 

Semi-major axis . 

a (au) 

0.04375 

±0.00070 ± 0.00023 

0.043 ± 0.002 

0.0435 ± 0.0006 

0.04473 ± 0.00024 

Inclination . 

i (degree) 

89.03 + 0.34 

89.8/°;! 

88.6/“;! 

89.138/!;®°! 

RV curve semi-amplitude . 

A'a (ms -1 ) 

82.7 + 4.2 “ 

82.1 ± 7.4 

74.5 ± 1.8 

76.16/!;!! 

Barycentric RV . 

7 (kms -1 ) 

4.8951 ±0.0021 

4.9077 + 0.0015 

3.95 ± 0.43 


Proj. spin-orbit angle. 

A (degree) 

7 + 5 





Notes. Where two errorbars are given, the first refers to the statistical uncertainties and the second to the systemat ic errors. 
a This value of K\ was determined from our RV data, which were taken during transit time only. The value reported bv lBakos et all 
(l2009f) . which is based on RV data with a much better coverage of the orbital phase, was therefore preferred for the determination 
of the other physical parameters of the system (see Sect0. 



a/R * 

Fig. 13. The projected orbital obliquity as a function of 
scaled o rbital distance, a/R+. for those systems with T e g < 
6150 K (lAnderson et al.l[2015l) . The empty circles represent 
near-circular orbits (e < 0.1 or consistent with zero), while 
the light-red circles depict eccentric orbits. The size of each 
circle is proportional to the corresponding pl anetary radius . 
The data have been taken from TEPcat (|Southworthll201 ih . 
The black points indicate the two TEP systems examined in 
this work, HAT-P-36 and WASP-ll/HAT-P-10. The error 
bars have been suppressed for clarity. 


time variation O-C of about 20 seconds in ten years that 
would be easy to measure with a space-borne photometer. 
On the othe r hand , if we assume Q* ~ 10 7 — 10 8 , as 
suggested bv lOgilvie fc Linl ( 2007 '). the orbital decay and 


alignment timescales are of the order of a few hundred 
Myr or a few Gyr, respectively. They are shorter than or 
comparable to the estimated age of the star. In conclusion, 
in this scenario the present alignment could result from 
the tidal interaction during the main-sequence lifetime, 
while the remaining lifetime of the system is at least of a 
few hundred Myr. In any case, HAT-P-36 is an interesting 
system to study the orbital and stellar spin evolution, as 
resulting from the tidal interaction and the magnetic wind 
braking of the star. 

For WASP-11, the relatively large semi-major axis and 
the small mass of the planet imply, on the other hand, a 
weak tidal interaction with an angular momentum exchange 
timescale that is comparable or longer than the age of the 
star. Therefore, the observed projected alignment of the 
system is likely to be of primordial origin. 


7. Summary and conclusions 

In the framework of the GAPS programme, we are observ¬ 
ing a sample of TEP systems hosting close-in hot Jupiters, 
with the HARPS-N spectrograph. The goal is to measure 
their spin-orbit alignment (by observing the RM effect) 
and study if this quantity is correlated with other physical 
properties. Here, we have presented an exhaustive study of 
the HAT-P-36 and WASP-11/HAT-P10 planetary systems, 
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Fig. 12. Sky-projected orbital obliquity as a function of several stellar, orbital and planetary parameters. Orange (white) 
circle s indicate systems in which the parent stars has an effective temperature higher (lower) than 6100 K (iDawsonl 
I2014H . The size of each cir cle is proportional to the corresponding planetary radius. The data have been taken from 
TEPcat (jSouthwortbl 1201 lh . The black points indicate the two TEP systems examined in this work, HAT-P-36 and 
WASP-ll/HAT-P-10. The error bars have been suppressed for clarity. 


based on data collected with four telescopes during transit 
events. 

By analysing these new photometric and spectroscopic 
data, we revise the physical parameters of these two sys¬ 
tems, finding that previous determinations agree with our 
more accurate results, within the uncertainties. 

Interestingly, we observed anomalies in three photomet¬ 
ric transit light curves of HAT-P-36 that, after appropri¬ 
ate modelling, turned out to be compatible with starspot 
complexes on the photosphere of the star. The main char¬ 
acteristics of the starspots were estimated and compared 
with those observed in other TEP systems. The HAT-P- 
36 starspot activity is confirmed by the analysis of the 
HARPS-N spectra, which give a stellar activity index of 
log Ajj K = —4.65dex, and by study of the whole photomet¬ 
ric data set collected by the HATNet survey. The frequency 
analysis of this time-series data revealed a clear modulation 
in the light curve caused by starspots, allowing us to get 
an accurate measurement of the rotational period of the 
star, P rot = 15.3 ± 0.4 days. A similar study, performed for 
WASP-ll/HAT-P-10 on both the WASP and HATNet sur¬ 


vey light curves, did not highlight any clear photometric 
modulation. 

The RM effect was partially covered in HAT-P-36 (be¬ 
cause of bad weather conditions) and wholly observed 
in WASP-ll/HAT-P-10. Thanks to the high-resolution 
HARPS-N spectra, the sky projection of the orbital obliq¬ 
uity was successfully measured for both the systems, indi¬ 
cating a good spin-orbit alignment. In particular, for the 
HAT-P-36 system, we were able to also estimate its real 
obliquity, obtaining ip = 25jz|i degrees. Our results are thus 
in agreement with the idea that stars with relatively cool 
photospheres should have small spin-orbit misalignment an¬ 
gles. However, looking at the entire sample of known TEP 
systems, for which we have an accurate estimation of A, 
there are several exceptions to thi s scenario. In this con¬ 
text, the case of HAT-P-18b (iEsposito et al.l[2014j) is em¬ 
blematic. The collection of more data is thus mandatory for 
disentangling the issue and understanding whether orbital 
obliquity really holds imprints of past migration processes 
that dramatically affected the evolution of giant planets. 

Acknowledgements. The HARPS-N instrument has been built by 
the HARPS-N Consortium, a collaboration between the Geneva 
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Fig. 14. True orbital obliquity as a function of the effective stellar temperature for 12 TEP systems. Data taken from 
TEPCat. The red point indicates the location of HAT-P-36 (this work). The temperature range is the same as that in 
the uppe r-left panel of Fig. HU The green-dotted line separates the cool TEP systems from those with T e g > 6100 K 
(lDawsonll2014 h 
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Appendix A: Supplementary tables 

The tables in this Appendix contain the detailed results 
of the predictions of different sets of stellar evolutionary 
models for the HAT-P-36 and WASP-ll/HAT-P-10 plan¬ 
etary systems. The final values of the physical parameters 
for each of the two systems (see Tables[S] and [9]), are calcu¬ 
lated by taking the unweighted mean of the five estimates 
of the different sets of model predictions (see SectO . 
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Table A.l. Derived physical properties of the HAT-P-36 system based on the prediction of different theoretical models. 



(Claret models) 

(Y 2 models) 

(BaSTI models) 

(VRSS models) 

(DSEP models) 

AT ( km s 1 ) 

196.19 ±4.95 

196.09 ±1.79 

193.00 ±1.62 

193.00 ±2.76 

196.17 ± 0.32 

M* (Mg) 

1.0503 ±0.0840 

1.0486 ±0.0287 

1.0000 ±0.0259 

1.0000 ±0.0426 

1.0500 ±0.0045 

R* (Re) 

1.048 ±0.029 

1.047 ±0.012 

1.031 ±0.013 

1.031 ±0.018 

1.048 ±0.010 

log P* (cgs) 

4.4192 ±0.0132 

4.4190 ±0.0100 

4.4121 ±0.0092 

4.4121 ±0.0109 

4.4192 ±0.0085 

M p ( JV/jup) 

1.877 ±0.130 

1.875 ±0.088 

1.816 ±0.087 

1.816 ±0.095 

1.876 ±0.079 

Rp ( ^Jup) 

1.313 ±0.038 

1.312 ±0.020 

1.291 ±0.021 

1.291 ±0.025 

1.313 ±0.018 

Pp (pjup) 

0.776 ±0.049 

0.776 ±0.046 

0.789 ±0.047 

0.789 ±0.048 

0.776 ±0.047 

e 

0.0654 ±0.0033 

0.0654 ±0.0030 

0.0665 ± 0.0030 

0.0665 ±0.0032 

0.0654 ±0.0030 

a (AU) 

0.024040 ± 0.000605 

0.024028 ±0.000219 

0.023651 ±0.000198 

0.023651 ±0.000338 

0.024038 ±0.000036 

Age (Gyr) 

4.3t?;t 

4 n +11 
+‘ U -0.8 

7 fi +4 ' 3 

5-0±|;g 

1 7+7-5 

A -* -0.1 


Notes. In each case g p = 27.0 ± 1.4ms 2 , p* = 0.913 ± 0.027 pg and T eq = 1788 ± 15K. 


Table A. 2 . Derived physical properties of the WASP-ll/HAT-P-10 system based on the prediction of different theo¬ 
retical models. 



(Claret models) 

(Y 2 models) 

(BaSTI models) 

(VRSS models) 

(DSEP models) 

Kh (km s 1 ) 

129.9 ±3.2 

127.7 ±2.0 

127.5 ±1.9 

127.5 ±2.3 

128.4 ±2.1 

M* (Mg) 

0.846 ±0.063 

0.804 ±0.039 

0.800 ±0.036 

0.800 ±0.043 

0.818 ±0.038 

R * (Re) 

0.785 ±0.023 

0.772 ±0.013 

0.771 ±0.015 

0.771 ±0.014 

0.777 ±0.012 

log 5 * (cgs) 

4.576 ±0.016 

4.569 ±0.016 

4.568 ±0.015 

4.568 ±0.018 

4.571 ±0.018 

M p (Mj up ) 

0.508 ±0.028 

0.492 ±0.020 

0.490 ±0.019 

0.490 ± 0.024 

0.497 ±0.023 

Rp ( -Rjup) 

1.006 ±0.031 

0.989 ±0.019 

0.987 ±0.020 

0.987 ±0.023 

0.995 ±0.022 

pp (pjup) 

0.467 ±0.028 

0.475 ±0.025 

0.476 ±0.026 

0.476 ±0.027 

0.472 ±0.026 

e 

0.0531 ±0.0021 

0.0540 ±0.0016 

0.0541 ±0.0017 

0.0541 ±0.0019 

0.0537 ±0.0019 

a (AU) 

0.04447 ±0.00111 

0.04373 ±0.00069 

0.04365 ±0.00065 

0.04365 ±0.00080 

0.04398 ±0.00070 

Age (Gyr) 

k 7+5.8 

1 -4.8 

8.9l 3 3 j 



6-2i ll 


Notes. In each case g p = 12.45 ± 0.50ms -2 , p* = 1.748 ± 0.074pg and T eq = 992 ± 14K. We were not able to constrain a 
reasonable value for the age of the WASP-ll/HAT-P-10 system by using the BaSTI and VRSS models, because the evolution of 
a star of the mass of WASP-11 is small for times shorter than the age of the Galaxy. 


18 






